Paraxial meridional ray tracing equations from the unified 
reflection-refraction law via geometric algebra 
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Abstract. We derive the paraxial meridional ray tracing 
equations from the unified reflection-refraction law using 
geometric algebra. This unified law states that the nor- 
mal vector to the interface is a rotation of the incident 
ray or of the refracted ray or of the reflected ray by an 
angle equal to the angle of incidence or of refraction. We 
obtain the finite meridional ray tracing equations by sim- 
ply equating the arguments of the exponential rotation 
operators. We then derive the paraxial limits of these 
equations with the help of sign function identities. We 
show that by embedding the sign functions in the ray 
tracing equations, we explicitly declare our chosen sign 
conventions in symbols and not in prose. 

1 Introduction 

In paraxial optics, it is customary to declare before- 
hand the adopted set of sign conventions, as done for 
example in Nussbaum and Philips [T]. But there are as 
many sign conventions as there matrix optics authors. 
And debate ensues [5] -[Sj. To get a taste of this contro- 
versy, let us quote Welford's 1974 critique of Conrady's 
convention [7] [8]: 

I have avoided the issue [of sign conventions] by 
simply using the universally accepted conventions 
of coordinate geometry, together with vectors and 
direction cosines, as a consistent system which 
agrees with what is done in other branches of 
physics. This conflicts with what has been taught 
at Imperial College on one small point: the sign 
of the paraxial convergence angle, u, has always 
been taken according to Conrady's convention but 
after due discussing with my colleagues I decided 
to reverse it so as to agree with the convention for 
direction cosines; the inconsistency could, we felt, 
no longer be justified. 



To resolve this sign problem, we propose the use of 
sign functions that take values of ±1, such as the three 
axial direction functions for the vector v: 

(1) 
(2) 
(3) 



These functions correspond to the signs of the direction 
cosines of a vector v. So by using these sign functions, 
we explicitly adopt Welford's Cartesian sign conventions. 

Another example of a sign function is the concavity 
function |9]: 
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where cr is the incident ray and r\ is the normal vector to 
the interface. If the interface is concave, then c^r; = 1; if 
convex, Co-j; = — 1. We used the concavity function before 
when we wrote down the unified reflection-refraction law 
in exponential form|10j: 



--(17] 



V 



(5) 

where cr' is the refracted ray and cr" is the reflected ray. 
From this unified law we derived the ray tracing equa- 
tions for finite and paraxial skew rays in spherical coordi- 
nates and for finite meridional rays in polar coordinates. 
The paraxial meridional rays we deemed then to require 
a separate treatment; we set it aside for a future work. 

In this paper, we shall continue our work. We shall 
start with a short review of geometric algebra and then 
proceed to geometric optics. We shall summarize the 
equations for finite skew and finite meridional ray trac- 
ing, and then use these to derive those for paraxial merid- 
ional rays in polar coordinates. We shall see how the use 
of sign functions makes the discussion of sign conventions 
unncessary. 
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2 Geometric Algebra 

In Clifford (geometric) algebra CZa.o the product of two 
vectors a and b is given by the Pauli identity [TT]: 



ab = a • b + i(a X b), 



(6) 



where i is the unit imaginary scalar. Note that the geo- 
metric product ab is an associative product, unlike the 
dot product a • b and the cross product a x b. 

The exponential function in geometric algebra is also 
well-defined: 



g 

cos|0| +i— sin|0 



(7) 



which is the generalization of Euler's theorem in complex 
analysis. If a is a vector perpendicular to 0, then we can 
show that 



a cos 101 — a X - — - sin \ 0\ 

' ' \e\ ' ' 



(8) 



Equation ([8]) states that ae'^ is the vector a rotated coun- 
terclockwise about the vector d by an angle |0|.[T2] 

3 Geometric Optics 
3.1 Finite Skew Rays 

Finite skew rays[T3] are the most general type of rays. 
The ray tracing equations for these rays are expressed in 
vector form . 

From its initial position Tq, a light particle travels by 
a distance s in the direction of the unit vector cr. The 
final position r of the light particle is 



r = ro + s(T. 



If at the position r, the outward normal unit vector to 
the interface is r] and the interface is spherical of radius 
R centered at C, then 



r = To + scr = C + Rt]. 



(10) 



After the light particle strikes the interface, the par- 
ticle is either be refracted or reflected. The directions 
cr' and <t" of the refracted and reflected vectors may be 
expressed in terms of the directions cr of the incident ray 
and (outward) normal vector r\ to the interface: 



cr 

a" 



(11) 

(12) 



where (3 and (3' are the angles of incidence and refraction. 



(3 — sin ^ |cr X r/j, 
13' — |cr X ryl, 



(13) 
(14) 



and Go-xr; is the rotational axis direction. 
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(15) 



We may also combine the laws of refraction and reflec- 
tion in Eqs. pT|) and (|12p into one, as given in Eq. ([5]). 
This unified law expresses the normal vector in terms 
of the rotations of the incident, refracted, and reflected 
rays about the vector e^xj) by angles f3, /?', and (3" — (3, 
respectively. The rotations are counterclockwise or clock- 
wise depending on the sign value of the concavity func- 
tion Cari- (Figure 1) 




Figure 1: The incident ray <t, refracted ray cr', 
and reflected ray cr". The rays make an angle of 
/3, /3', and /3" with respect to the normal vector rj, 
respectively. The interface is convex, Carj = — 1. 



(9) 3.2 Finite Meridional Rays 



Finite meridional rays |14j are rays that lie on the same 
plane. Here, we choose this plane to be the zx— plane, 
with the z— axis along 63 as the optical axis. 

Let us define the incident, refracted, reflected, and nor- 
mal vectors as vectors in za;-plane: 

(16) 
(17) 
(18) 
(19) 

where is a counterclockwise rotation angle measured 
from 83. 

If we also define the position vectors r and rp as vectors 
in the za;— plane, 
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(20) 
(21) 
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and the center 



C = Zees 



(22) 



of the interface to be along the optical axis 63, then 
Eq. (fTU]) separates into 

z — zo + scos6„ = C + Rcos9^, (23) 
X = xq + s sinOg^ — RsinOfj. (24) 

Equations (^5)1 and are the ray propagation equa- 
tions for finite meridional rays. (Figure 2) 

On the other hand, substituting the definitions in 
Eqs. (Uni) to ([TO]) back to the unified reflection-refraction 
law in Eq. ([5]) and employing the identities 



„ _ ie2(c„„-l)7r/2 



we arrive at 



2(c<T,, - l) + 6'r, 



(25) 
(26) 
(27) 



Equation ([28)1 is the unified reflection-refraction law for 
finite meridional rays. We can show that this equation 
contains a restatement of the Bessel-Conrady refraction 
invariant, by replacing 9^ by U and CcrnC(„xi-i)yP by I. 



s 



















Figure 2: A ray travels a distance s until it in- 
tersects an interfaces of radius R. 



3.3 Paraxial Meridional Rays 

Paraxial meridional rays are rays that lie on the same 
plane (zo;— plane) and make a small angle with respect 
to the optical axis 63. Mathematically, we say that if v 
is a paraxial meridional ray, then the polar angle 9^ of 



V may be expressed in terms of the small positive angle 
9vz that V makes with the optical axis [17]: 



-(1 C^^) ~l~ Cvz^vx9vz : 



(29) 



where 9^ is the polar angle of v. We can easily verify 
that 



cos 9.U — Cv 



(30) 
(31) 



Note that Cyx and Cyz are the axial direction functions 
defined in Eqs. ^ and 

3.3.1 Propagation 

Using the approximations in Eqs. ([50)1 and (|?T|) . Eqs. 
and (|24p simplify to 

Z = Zq + C^zS ^ C + Cr,zR, (32) 

X = xo + Ccrxs9^z ^ CjjxR9,,z- (33) 

Equations ([5^ and ([55]) are the position-height relations 
for paraxial meridional ray propagation. 

For the polar angle 9^, we know that it is conserved 
during propagation or translation: 



(34) 



Using Eq. (^5]) . Eq. ([M)) may be expanded as 



-^(i-Ccrz)+Ca-zCcrx9az = -^[i-CaQz)+CaOzCaOx9aOz, (35) 

Because the light particle moves in the same direction of 
the z-axis during propagation, then 



^(Tz — ^aOz 

so that Eq. l[35|l reduces to 

C(7x9(JZ (^(TQx9(j{)Z- 



(36) 
(37) 



Thus, the inclination angle of a ray from the optical axis 
remains invariant under ray propagation. 

Substituting Eq. (157)) back to the propagation equation 
in Eq. (|55)l . we obtain 



a; = Xo + SCcrQx9crOz- 



(38) 



Equation ()55)) expresses the height a; of a ray from the 
optical axis as a function of the ray's initial height xo, 
distance travelled s, and direction angle ^o-oz. Notice that 
the magnitude 9„qz of the paraxial angle is divorced from 
its x-direction sign function c„qx- 
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3.3.2 Refraction 

The refraction law may be extracted from the unified law 
in Eq. (|28ll : 



(39) 



We know that the angles of incidence and refraction are 
related by Descartes-Snell's law: 



n'sin/3' — nsm(3, 



(40) 



where n and n' are the refractive indices of the medium 
containing the incident and refracted rays, respectively. 
In the paraxial limit, Eq. (j40| becomes 



Hence, 



where 



n'f3' = n(3. 



n 



(41) 

(42) 
(43) 



Equation is the paraxial refraction law in terms of 
the polar angle 9^ of the incident ray and the angle of 
incidence f3. 

Using the unified refraction-reflection law in Eq. (|28p . 
we solve for the angle of incidence (3 in terms of the polar 
angle 9ri of the normal vector: 

7r 

CCTr,C(^xr,)a/3 = -(Cff^ - 1) + 6*^ - 6*^, (44) 

Substituting this back to Eq. and rearranging the 
terms, we arrive at 

Oa' = f^O, + (1 - M)(^(c.r, - 1) + 0^,). (45) 

In paraxial approximation, this is 

"^(1 ^{t' ^{t' z^a' x9a' z 

= yu(l - Ccr^) - + l-J.CazCax()crz 
+ (1 -M)(Cff>) 



(46) 



Now, we can verify that the following sign identities 
hold: 



Car] 



(47) 
(48) 



The first equation states that the relative directions of 
the incident and refracted rays with respect to the optical 
axis are the same; the second states that the concavity 
function is equal to the product of the relative directions 



of the incident ray and normal vector with respect to the 
optical axis. 

Employing the sign identities in Eqs. (|47|) and (|48p. 
Eq. ((46|l reduces to 



Ca'xOa'z = m' + flCax0az + (1 " y)c„r)Cy^xd-nz , (49) 

where 

7r 

m' = {I- ^J.){criz - c^zCnz - Caz + l)-^- (50) 

Let us analyze the angular function m' in Eq. (|50p by 
considering two cases: 

m' = 0; = +1, (51) 

m' = 7r(l-^)(l + c^^); = (52) 

Because the direction c^^ of the normal vector with re- 
spect to the 63 is arbitrary, then m' 7^ in general for 
backward propagating rays (caz = ~1)- This is an in- 
convenient case. So we shall leave this for a future work 
and impose that 



^az — ^a' z 



-1 



(53) 



for all our equations. Thus, Eq. (|49|) reduces further to 

Ccr'xdcr'z = ^J-CaxOaz + (1 ^ fJ-)Cr,zCr,xOrjz ■ (54) 

Equation (|54p is the refraction law for forward propagat- 
ing paraxial rays. 

Using the result in Eq. ([33]) . 

C-qx^^-qz — 75 J (55) 
it 



Eq. ([54]) becomes 



Ca'xda'z = tJ'CaxOaz + (1 - /i)c,,z — . (56) 

-ft 



Multiplying this by n' yields 



where 



^ct' x^ ct' z ^C(T.r^rT2 Px. 



P = Cr, 



(57) 



(58) 



R 

is the power of the interface. Notice that the radius R 
of the interface is always positive; the z— axis direction 
function Cr,z of the normal vector tj takes care of the sign 
traditionally possessed by R. 

3.3.3 Reflection 

From the unified refraction-refiection law in Eq. (I28p . we 
get 

9a" = -n + 6a + 2Ccrr,C(crxj))2,/3- (59) 
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Using the result in Eq. (|44| . Eq. (|59|) becomes 

e^" = (c^„ - 2)^ - 0, + 26.,. (60) 
In the paraxial limit, Eq. (|60|l reduces to 

TT 

■^(l-Co-''^) + Ccr"zCCT":r6'o-"z (61) 

TT 

= 7r(c^^ - - 1) + -(c^^ - 1) 

^az^ax^az ~l" '^^'qz^rix^'qz- (^2) 

Because the z— direction of the reflected ray is opposite 
to that of the incident ray, then 

Ca"z = -Caz, (63) 

SO that Eq. becomes 

^ ^ax^az ^" '^^trrj^rix^riz : (64) 

after multiplying by c^^ and using the sign identity in 
Eq. jUl). 

Let us analyze the angular function 

m" = 7r(-c,,2 + c^zCrj^ + 2c„z) (65) 

by considering two cases: 

to" = 27r; c,^^ = +1, (66) 

to" = -27r(c,,, + 1) = {0,-47r}; c,, = -1,(67) 

since — ±1. Because 47r = 27r = 0, then we may 
simply set m" = 0, so that Eq. ([64]) simplifies further to 

^a" x^a" z — ^(Tx^az '^C{7r]^r]x^r]z ■ (68) 

Employing the identity in Eq. ([55]) . Eq. (|68p becomes 

^(j"x^a"z — ^ax^az 2Ccrj^— . (69) 

il 

If we impose that the incident ray moves along 63 , then 
Caz = +1, so that Eq. ([55)1 simphfies to 

X 

^a" x^cr" z ^ax^crz ^^rjz~7^- ('^^) 

H 

Multiplying Eq. ([TDD by n" = n, 

n"Ca"xda"z = nCaxOaz " 2nc^2^, (71) 

H 

and comparing the result with the refraction law in 
Eq. (|57p . we see that we may express the mirror power 
P" as 

^'" = c,.f. (72) 

Notice that except for the factor of 2, the power P of a 
mirror Eq. (172^ is similar to that of a lens in Eq. (I58p . 
which is what we expect. 



4 Conclusions 

In this paper, we used geometric algebra to derive from 
the unified reflection-refraction law in exponential form 
the paraxial meridional ray tracing equations in polar 
form, by equating the arguments of the exponentials and 
employing the properties of sign functions. The sign 
functions, such as the concavity and axial direction func- 
tions, make explicit the Cartesian sign convention used, 
though in symbols and not in words. We hope that these 
sign functions would be universally adopted to finally 
settle the never-ending debate on sign conventions. 
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